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We perform lattice calculations of the spatial 't Hooft ^-string tensions, ojt, in the deconfined 
phase of SU(N) gauge theories for N = 2,3,4,6. These equal (up to a factor of T) the surface 
tensions of the domain walls between the corresponding (Euclidean) deconfined phases. For 
T ^> T c our results match on to the known perturbative result, which exhibits Casimir Scaling, 
Ojt x k{N —k). At lower T the coupling becomes stronger and, not surprisingly, our calculations 
show large deviations from the perturbative T-dependence. Despite this we find that the behaviour 
d&k/dT oc k(N — k) persists very accurately down to temperatures very close to T c . Thus the 
Casimir Scaling of the 't Hooft tension appears to be a 'universal' feature that is more general 
than its origin in the low order high-T perturbative calculation. We observe the 'wetting' of these 
fc-walls at T ~ T c . Our calculations show that as T — * T c the magnitude of &k(T) decreases rapidly. 
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1. Introduction 

In the Euclidean formulation of SU(A0 gauge theories at finite T , deconfinement is associated 
with the spontaneous breaking of a Zn centre symmetry. The obvious order parameter is the trace 
of the Polyakov loop, l p . In the deconfined phase the effective potential for the Polyakov loop trace 
averaged over the volume, l p , has its minimum at l p ^0, so the symmetry is spontaneously broken 
and there are N possible deconfined phases. If f p « 3 = e l7ll N we label the phase by k. When 
two of these phases, k\ and ki say, co-exist they will be separated by a domain wall whose surface 
tension will depend on k = k\ — &2, as well as on ,/V and T. The tension of these domain walls 
is equal, up to a factor of T, to the spatial 't Hooft string tension, which is in principle related to 
confinement. 

At high T one can calculate in perturbation theory. To two loops one finds ^ |3j] : 



where C2 — 0.09. 

The factor k{N — k) is the k dependence of the Casimir, Tr z ^T a T a , where M is the totally 
antisymmetric representation of a product of k fundamentals of SU(A0. It is the factor one obtains 
when calculating the Coulomb interaction between sources in such a representation. In D=l+1 
SU(A0 gauge theories, the tension of the confining ^-string that connects such sources has precisely 
this dependence. There are old speculations [Q] that this 'Casimir Scaling' holds in D=3+l and 
numerical calculations show that it is a good (but not exact) approximation & A question we will 
try to answer in this paper is whether the Casimir Scaling in eqn( |l.l[ ) survives at very much lower 
T . This question connects to the role of Casimir Scaling in confinement, since these domain walls 
are closely related to the centre vortices that provide a possible mechanism for confinement!! E0- 

2. Preliminaries 

We discretise Euclidean space-time to a periodic hypercubic lattice. The jit = direction is the 
temperature direction and the domain wall spans the L\ x L2 torus. 

Our order parameter will be based on the Polyakov loop l p . Above T c there are N degenerate 
phases in which (l p ) = Zkc(fi) where c(j8) is a real- valued renormalisation factor. They can co-exist 
at any T >T C and will be separated by domain walls. If l p in the two phases differs by a factor Zk 
we refer to the domain wall as a &-wall. 

To study a &-wall we use a 'twisted' plaquette action to enforce the presence of a single domain 
wall. The twisted action is defined by 



o^ = k(N-k) 



An 2 T 3 



{l-c 2 g 2 (T)N} 



3V3 VgHf)N 



(1.1) 




(2.1) 



where z(p) = 1 for all plaquettes except 



z(p = {nv,x ll })=z k = e 



jUV = 03; *o = xo, *3 =X3, x h2 = l,...,L li2 . 



(2.2) 
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That is to say, the plaquettes in the entire (0,3)-plane at xo = xq' and X3 = x^ are multiplied by 
Zk £ Zjy. The Polyakov loops on either side of the plane will differ by a factor of Zk- Periodicity 
in X3 then demands that at some X3 the Polyakov loop must suffer a compensating factor of Zk^ - a 
domain wall. Thus we ensure that each configuration possesses at least one fc-wall. 

We calculate the average action with and without a fc-twist as defined above. The difference is 

ac-/c\ /c \ _ dlnZ o dlnZfc _ d F k -F _ d a^A 

AS k = (S k ) - (S ) _ - - (2-3) 

where is the surface tension of the domain wall and A = a 2 L\L2 is its area. 
We see from eqn(|l.lD that we expect 

d o^A T ^ MN-k) A% 2 UL 2 d 1 

when leading order perturbation theory is accurate. The factor of a(Lo) contains the 0(a 2 T 2 ) = 
0(1 /Lq) lattice correction, obtained by a numerical evaluation of the expression given in jj]. 

3. Results 

3.1 Surface tension 

Our most extensive results are for SU(4) with Lo = 4. We show the values of AS^T 1 for 
these calculations in Fig. [T], compared to the two-loop perturbative expectations using a mean-field 
improved coupling and the bare coupling. We also include a 'good' coupling evaluated at low 
energy scales: the Schrodinger functional coupling which has been calculated non-perturbatively 



in [JTTJ]. We plot the ratios AS^/ASj, in Fig. |. 

We obtained results at one additional point, not included in the graphs, at j3 = 20, correspond- 
ing to T ~ 1000r c . Here we found excellent agreement with perturbation theory, as expected. 

In Table | we list our results for AS*, for SU(4) with L = 5 and for SU(2), SU(3), and SU(6). 



We include our results for SU(4) with Lo = 4 at the same temperatures for comparison 

'IV 



For SU(4) with Lq = 4, as T — > T c we find AS^ grows much more rapidly than the perturbative 



prediction, reaching a factor of 18 at T ~ 1.02r c . This tells us that da^/dT becomes much 
larger than the low-order perturbative expectation as T —> T c , implying that o^(T) is increasingly 
suppressed relative to its perturbative value as we approach T c . The ratio of the AS 1 ^ satisfies 
Casimir Scaling to good accuracy, though there is some evidence for a discrepancy below 1.5T C . 

To investigate the continuum limit we compare to L = 5. The discrepancy with perturbation 
theory is the same at T ~ 1.88T C and similar at T ~ 1.02T C . It is clear that a large and growing 
mismatch with perturbation theory as T — > T c is a feature of the continuum theory. The ratio of the 
ASw continues to be close to Casimir Scaling, so that is also a property of the continuum theory. 

For SU(6) we observe at both values of T precisely the same discrepancy with perturbation 
theory as we saw for SU(4) at the same value of Lo. In addition the AS{y ratios continue to satisfy 
Casimir Scaling. Taken together this tells us that the derivative of the domain wall tension has 
no factors of k and Af except for the Casimir scaling factor k(N — k) and its dependence on the 
't Hooft coupling, g 2 N. Our SU(2) and SU(3) calculations show a very similar discrepancy with 
perturbation theory. Thus the suppression of o^(T) as T — > T c is largely independent of N. 
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To estimate the suppression of a^, we can in principle calculate it by interpolating 
in p and then integrating from large /3 down to the desired value of T, but our calculations are not 
dense enough in /3 for this. We can obtain a qualitative picture by assuming some functional form 
for do^/dp and fitting it to the calculated values. Choosing the L = 4 SU(4) calculation, we 
make a fit using a simple modification of the one-loop formula for the /c-wall free energy: 

= [1 loop] +aexp(-ZV/3/-/3/ f ) (3.1) 

Using this fit we obtain the T-dependence of Gy, shown in Fig. ^[ We display the uncertainty from 
the errors in the fitted parameters, but it is clear that we cannot reliably estimate the systematic error 
inherent in the choice of fitting function. The qualitative picture is that is strongly suppressed 
at T < l.5T c . This suggests that it may reach zero at some 7# somewhat below T c , causing a 
second-order phase transition due to the condensation of spatial 'tHooft strings. 

3.2 Profile 

We can average the Polyakov loop over the transverse coordinates to obtain the profile l p (x^). 
Comparing with the perturbative prediction rf2L ||] we find good agreement down to T ~ 1.88r c . 
However, at T ~ 1 .02T C the profile is very far from the perturbative expectation. This observation 
makes it all the more remarkable that we continue to see Casimir Scaling at such very low T. 

3.3 Wetting 

A fc-wall can interpolate between two deconfined phases by passing through the origin of the 
complex l p (n) plane, i.e. the confined phase. Whether this will happen or not depends on the 
relevant surface tensions. We investigated this by a series of runs very near T c in SU(4) and SU(6). 
We found that all the walls split into a pair of connned-deconfmed walls ('wetting') over a range 
of ~ 0.0ir c in T. This small range occurs because the latent heat is much larger than the domain 
wall tensions, which suppresses the breaking up unless T is extremely close to T c . 

4. Conclusions 

We have shown that surface tensions of domain walls seperating different deconfined phases 
are close to Casimir Scaling, i.e. 

ok=k(N-k)f(g 2 (T)N,T/T c )T 3 , (4.1) 

to a good approximation. The surface tensions are strongly suppressed as T — > T c . The domain wall 
profile is far from the perturbative prediction near T c . We also observed 'wetting' of the domain 
walls very near to T c . 
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Table 1: Results for AS* . The prediction is the two-loop perturbative expectation, using a mean-field 
improved coupling. 
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Figure 1: Action per unit area of the k = 1 domain wall in SU(4) with aT = 0.25. Monte Carlo values, +, 
compared with perturbation theory based on various couplings: g 2 (a), solid line, gj(a), long dashed line, 
gs F (T), short dashed line. 




Figure 2: Surface tension in units of T, using the functional form in eqn(3.1) fitted to our results. For 
comparison we show the 1-loop (solid line) and 2-loop (dot-dashed line) perturbative results using a mean 
field improved coupling. All for the k = 1 wall in SU(4). 
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Figure 3: Ratio ASl/AS l w in SU(4). 



179/6 



